Neutral Z boson pair production due to radion resonance in the 
Randall- Sundrum model: prospects at the CERN LHC. 

Prasanta Kumar Das ^ 

Chung- Yuan Christian University, 
Chung-Li, Taiwan 320, Repubhc Of China 

Abstract 

The Neutral Z boson pair production due to radion resonance at the Large Hadron Collider 
(LHC) is an interesting process to explore the notion of warped geometry (Randall- Sundrum 
model). Because of the enhanced coupling of radion with a pair of gluons due to trace enomaly 
and top(quark) loop, the radion can provide larger event rate possibility as compared to any 
New Physics effect. Using the proper radion-top-antitop (with the quarks being off-shell) 
coupling, we obtain the correct radion production rate at LHC and explore several features of 
a heavier radion decaying into a pair of real Z bosons which subsequently decays into charged 
4 l{l = e, fi) leptons (the gold-plated mode). Using the signal and background event rate, we 
obtain bounds on radion mass and radion vev (0) at the 5a, 10a discovery level. 
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1 Introduction 



The Randall-Sundrum (RS) model [1], proposed as a resolution of the well-known electroweak 
hierarchy problem, is particularly interesting from the phenomenological point of view [2]. 
According to this model, the world is 5-dimensional and the extra spatial dimension is jZ^ 
orbifold. The metric describing such a world can be written as 

ds^ = e-^'^'^^^^^^cix'^dx" - Rlde\ (1) 

where k is the bulk curvature constant and R^. determines the size of the extra dimension. 
The variable 9 parametrizes extra dimension. The model is constructed out of two branes 
which arc located at two orbifold fixed points = and 9 = n, respectively. The brane 
located at = (where gravity peaks) is known as the Planck brane, while the brane located 
at ^ = TT [where the Standard Model (SM) fields reside and gravity is weak] is known as the 
TeV brane. The factor e~^*^=l^l appearing in the metric is known as the warp factor. The 
length Rc, the distance between the two branes, can be related to the vacuum expectation 
value (vev) of some modulus field T{x) which corresponds to the fiuctuations of the metric 
over the background geometry given by Rc- Replacing Rc by T{x), we can rewrite the RS 
metric at the orbifold point ^ = tt as 

^ gjiJdx^dx" (2) 

where g'^J — e~'^^'''^^^^r)ni, — (^^) Vni^- Here p = and M5 is the 5-dimensional Planck 
scale. One is thus left with a scalar field (j){x) = — (0)] which is known as the 

radion field [3] . Golberger and Wise [4] proposed a mechanism for generating the potential 
for the field T{x), with its minima at Rc which satisfies kRc ~ 11 ~ 12, a requirement for 
the hierarchy resolution. With this nonminimal version of the RS model with the radion 
0(a;) stabilzed at (0) via the Golberger and Wise mechanism, one can do a lot of interesting 
phenomenological studies. In particular, the radion, which can be lighter than the other 
low-lying gravitonic degrees of freedom in the RS model, will reveal itself first directly in 
the collider experiment or indirectly in the precision measurement and verify our notion of 
extra dimensions. Studies based on observable consequences of radion are available in the 
literature [5], [6]. 
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Here we want to look at the pair production of Z bosons mediated by the radion at large 
hadron collider (LHC). The crucial roles in enhancing the signal in this channel are (a) the 
accessibility of the radion resonance for > 2mz and (b) the relatively enhanced radion 
coupling (due to trace anomaly and top quark loop) with a pair of gluons at LHC energies, 
which we will discuss in succession in later sections. Finding the radion via the above channel, 
the so-called gold-plated mode, is already available in the literatures. However some error has 
been translated in those studies due to the improper treatment of the Feynman rule in the 
radion coupling with a pair of off-shell top quarks. In the earlier studies, while estimating 
the radion production rate due to gluon fusion (i.e. in LHC context), radion coupling to a 
pair of top quarks is assumed to be proportinal to the mass rrit of the top quark, which is 
definitely true if the produced top quarks are on-shell and obviously does not hold in the 
present case, since the top quarks are appearing inside the loop i.e. they are off-shell. This 
is the main difference between us and the already existed literatures and will be discussed 
in detail in Section 2. 

The organization of the paper is as follows. In Sec. 2, we discuss the various interactions 
of the radion with the SM fields and determine the effective interaction of radion {(f)) with a 
pair of gluons [g) which is due to the top loop and QCD trace anomaly. We discussed the 
difference in the production rate obtained by us with those alreay existed in the literatures. 
The general features of Z boson pair production via radion are discussed in Sec. 3. Sec. 
4 contain discussions of the predicted signals for m<^ > 2mz, which is followed by the 
background estimate. We conclude in Sec. 5 



2 Effective interaction of radion with the SM fields 

Radion interactions with the SM fields confined on the TeV branc are governed by 4- 
dimensional general covariancc. It couples to the trace of the energy-momentum tensor 
of the SM fields in the following manner: 

^^nt = -^fl^iSM), (3) 
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where (0) is the radion vev. The trace of the energy-momentum tensor of the SM fields is 
given by 



+{2mlh^ - di^hd^h) + ■■■ (4) 



The photon and the gluons couple to the radion via the usual top loop diagrams. It is 
important to note at this point that as long as the fermions (produced from the decay of 
radion) are on-shell, their couplings are exactly same as the SM Higgs boson i.e. proportinal 
to the fermion mass and this makes it very difficult to isolate radion from the Higgs boson by 
studying the final state decay products. However, if they are off-shell (e.g. radion production 
via top loop in gluons fusion), the situation is different due to the kinetic term of (j) — t — t 
coupling term in the kinetic energy part of the Lagrangian (see the discussion in section 2.1). 
Considering the correct Feynman rule for (j) — t — t coupling [ i.e. the kinetic term (missing 
in the literature) along with the mass term (present in the literature) ] and thus finding 
the correct radion production due to top loop in gluon fusion in the present conteoct, differs 
crucially from those already existed in the literature [5] . Besides the top loop, an added source 
of enhancement of radion production in gluon fushion is the QCD trace anomaly, which we 
discuss below. For the sake of completeness, we first derive the lagrangian comprising radion- 
top coupling and then discuss the effective radion-gluon coupling which is due to top loop 
and trace anomaly. 



2.1 Lagrangian for the radion-top coupling 

The radion coupling to the top quark in the Randall- Sundrum model can be derived from 
the following action 

(5) 

Note that the action (Eq. (5)) should also contain the Yukawa term -^Hi/ji/j with H being 
the SM Higgs field. Such a Yukawa term is irrelevant for the present study and we will not 
consider this any further. The field e"'^ is the contravariant vierbein field for the visible brane 
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and g^,is is the determinant of the metric g^]^ (Eq. (2)). In the presence of radion fluctuation 
it satifies the normalization condition 

(6) 



is the covariant derivative with respect to general coordinate transformation and is given 

by 

is the spin connection and it can be computed from the vierbein fields by using the 
relation 



where Y^ah is given by the expression = | [7a, Tfc]- In the presence of radion fluctuations 
on the visible brane, the spin connection is given by 



w 



ah 



(8) 



The covariant derivative of the fermion field then becomes 

where the 7^ are independent of space time coordinates. The action that determines the 
radion couplings to top quark can therefore be written as 



5i 



/ 



J. 



I 



(0). 



= / d'^x 



+ J d^x 



+ 



(9) 



We keep terms containing one field [ the second term in the last line of Eq. (9)] which 
gives rise to the radion-top interaction with the following Feynman rule: 
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FIG. 1. Feynman diagrams for the radion coupling with a pair of top quarks. 



i [3 



(10) 



(0) L2 



where p, correspond to four-momenta of t and t quarks (see Fig. 1) and rrit, the top 
quark mass. Note that if the top quarks are on-sheU, the above Feynman rule can be shown 
(using the equation of motion) to be equal to similar to the corresponding SM liiggs 

boson coupling: —i^ where the electro weak vev v = 247 GeV. These two couplings are 
exactly identical in the limit (0) —>■ v. In Eq. (9) the Ellipsis (• • •) corresponds to terms 
involving two,three fields [see [7] for the details]. Note that in the above ijj — (X) ^ ■0, 



and rrit — (-^j "^t- In the following, we shall assume that all fields and parameters have 
been properly scaled so as to corresponds to the TeV scale and drop the tilde sign. 

2.2 Effective radion-gluon coupling: QCD anomaly and top loop 
2.2.1 QCD anomaly at order ag 

The radion couples to the trace of the energy- momentum tensor (T^jy), which can be equated 
with the four-divergence of the dilatation current (I?^) associated with scale symmetry of 
the theory i.e. 



Scale symmetry is a good symmetry at the classical level (i.e. d^D^ = = 0) in the limit of 
zero mass and in the absence of dimensionful couplings. However, quantum correction breaks 
the scale invariance and the divergence of the dilatation current receives the anomalous 
contribution given by d'^D^ — Tj^ ^ 0. Such a non-zero anomalous generated due to 
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the quantum correction (zero at the classical level), is known as the Trace Anom,aly or Weyl 
Anomaly term and it can be shown to be proportional to the relevant beta function of the 
underlying gauge theory [8]. The radion coupling to the massless gluons (or photons) is 
crucially controlled by such a trace anomaly term. For gluons, the interaction lagrangian 
due to this trace anomaly term can be written as 

r^(5M)~ = ^ ^^^G^G"^'' , (11) 

where Qs is the strong coupling constant and (3s{gs)/2gs — —[as/Sir] hqco- Here hqcD — 
11 — 2n//3 and is the number of quark flavours. 

2.2.2 Top loop at order ag 

At order ccg, a radion can couple to a pair of gluons via the one- loop diagrams of Fig. 2. 
We follow the usual [9] procedure to evaluate the amphtude for the process — > gg. Consis- 
tency with the Lorentz invariance and the transverse nature of the gluon [i.e. P'^e'^*{P) — 
Q^e^*{Q) = 0] leads to the following form of the amplitude: 

F^{e)P''Q^ + F2{e)r)^''] £^*(P)£^*(Q) (12) 

with rj^i, = diag{l, ^1, —1, and k = P -\- Q. Imposition of the current conservation 
Qi.M^'' = (follows from the gauge symmetry) relates the two form factors as 

F2{e) = -p.Q F,{e). (13) 

Inserting this in Eq.(12), we find the amlitude A1(0 — > gg) as 

M(0 ^ gg) = F^(k') [P^Q^ - P-Qfj^ s;*(P)sl*(Q). (14) 
We now evaluate the set of diagrams of Fig. 2 and find Fi{k^) as 

-Iqcd, (15) 



as{fI^)SabQ} 



27r(0) 

where Q/ is 2/3 for the top quark. Diagrams 2(c) and 2(d) do not give rise to contributions 
of the form Eq. (14), a necessary requirement which follows from the gauge invariance and 
other constraints (discussed above). Thus their contributions to the effective (f) — g — g 
interaction, are simply zero. Only diagrams 2(a) and 2(b) contribute. 







FIG. 2. Feynman diagrams with a top quark loop contributing to the process 
The loop function Iqcd arising in the process is given by 



gg at order as 



iQCD 



20 - 12x/^7^ X{xt) - 2{Axt - l)\\xt)\ for Xt > 1, 

20 - UVT^t K^t) - 2{1 - ixt)X\xt)] forxt<l, (16) 



where Xt = Amf/m'^ and is the radion mass. Note the presence of the constant term 
and the term hnear in X{xt) arising in the loop function. They, found to be absent in the 
literature [5] , arises due to the presence of the kinetic term in the (f) — t — t coupling. As we 
will see, this results in an important role in the radion production rate. The function \{xt) 
is given by 



A(xt) 



tan 



In 



^ ^Xt - I) 



for Xt > 1, 

for Xt < 1. 



(17) 



l-^/^^^ 

Note that for Xt > 1, n/ = 5, so 6qcd [arising in the trace anomaly term (see Sec. 2)] is 
23/3, while for Xt < 1, nf — 6 and Bqcd — 21/3. Note that the Xt < 1 i.e. > 2mt, 
corresponds to the top resonance inside the loop and is the region where the heavy quark 
appriximation breaks down. 
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On the whole, the effective (p—g—g interaction finds contributions from the trace anomaly 
and the top quark loop and hence, can be written as 



M<t>) 



bQCD + QpQCD [P^'Q^-P.QV 



(18) 



With the above interactions, we calculate the decay width r[0 — > gg] as 



99] 



JQCD 



f^QCDl ■ 



327r3((/))2 

Radion decay to other SM fields have been discussed in Refs. [5, 10]. 



(19) 



3 Neutral Z boson pair production: LHC prospects 

The neutral Z boson pair production from gluon-gluon fusion at the LHC has been exten- 
sively studied in the context of the SM as well as several of its extensions such as the minimal 
supersymmetric standard model (MSSM). A precise estimate for the Z boson pair production 
due to Higgs resonance and a detailed analysis of it's subsequent decay into charged leptonic 
final states is available in the hterature [11]. Any excess over this can be interpreted as the 
signature of new physics [12]. The mediation of the heavier neutral Higgs in supersymmetry 
has been shown to be the source of some enhancement in a region of the parameter space. 
Theories with extra dimensions, too, have been explored in this context, both in the ADD 
(Arkani-Hamed, Dimopoulos and Dvali [13]) and RS models [14]. It has been reported in 
these papers that the mediation of gravitons can boost the Z boson pair production rates. In 
this paper we wish to emphasize that the presence of a stablized radion in the nonminimal 
RS context is particularly significant. This is because (i) whereas the graviton resonance 
in the RS scenario is usually at too high a mass to be significant, a relatively less massive 
radion can be within the kinematic reach, and also (ii) the enhancement of radion coupling 
to a pair of gluons via the trace anomaly term and top loop picks up the contributions. We 
have computed the predicted rates for the case > 2mz and analyzed the viability of the 
resulting signals with appropriate event selection strategics. 

Using Breit-Wigncr approximation for the resonant production of a radion, one finds the fol- 
lowing expression for the cross section for pp ^ — > ZZ (the dominant partonic subprocess 
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being gg ^ (f) ^ ZZ): 

(7zzipp ^ ^ ZZ) = j dxi J dx2 gi{xi, Q^) g2{x2, Q'^)^zz{9g ^ (j) ^ ZZ), (20) 



where 



(Tzzygg ^ (p ^ = t^t-t^ — ^7-77 Pqcd + V/Jqcd -ttt jTi"^ — 27^ 

2048 7r'^(0)* •' [(s - m^)^ + m^r^J 

ri'/'^fi'fi'] s[s^ — 4m|s + 12m|] 



64 (0)2 p - m2)2 + rri^r^] ' ^^^^ 
In obtaining Eq.(21), we made use the decay width of radion (0) into a pair of gluons {g) 
i.e. r[0 — > gg] and assumes the vahdity of time reversal symmetry. Here s {— xiX2s) is the 
centre-of-mass energy for the partonic subprocess, ^/s{— 14 TeV) corresponds to the proton- 
proton centre-of-mass energy, gi{xi,Q^) and g2{x2,Q^) are the gluon distribution functions 
of the two coUiding protons, xi and X2 are the momentum fractions of the partons (gluons) 
of the colliding protons and m<^ and denotes the mass and total decay width of the radion. 
The factorization scale and renormalization scale //^ appearing in g{x,Q^) and as(//^) 
are fixed at = /i^ = s in our analysis. We have used CTEQ4L [15] parton distribution 
functions, setting the renormalisation scale at the partonic subprocess centrc-of-mass energy 
(VS)- We have checked that the predicted results are more or less unchanged if this scale is 
set at radion mass (m^). 

In Eq. (21), we find that away from the radion resonance, the term proportional to is 
of little consequence, and the total rate falls as 1/(0)^. This is because the decay width of 
the radion in any channel is proportional to 1/(0)^. Near resonance, on the other hand, the 
contribution is dictated by the term with in the radion propagator. Since, in the narrow 
width approximation, one can write l/[(s — m^)^ + m^F^] ~ tt6{s — m'^)/m^T^, the overall 
rate falls as 1/(0)^ near resonance. 



4 Numerical Analysis 

4.1 Event selection and results: > 2'mz 

The signal of Z-pair production due to radion resonance depends on the final states produced 
by the decay of Z-boson. In our analysis, we consider the mass range > 2 mz and in 
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such a case, the BR^cj) — > bb) faUs to 0.0028 i.e. ~ 3% (say at = 200 GeV) and hence the 
bb mode for detecting the radion resonance, is not worthwhile to look. On the other hand for 
the same value, the BR{(f) — > ZZ) becomes 13%, and can be considered as quite useful 
for detecting the radion resonance. We will consider the scenerio in which the produced Z 
will decay mainly leptonicaly (i.e. I — e,ii). So, the final state contains 4 charged leptons 
and normally such events may have appreciable SM background arising from the following 
sources: 

• The electro weak production of — e,ij), 

• Higgs production due to gluons fusion, which after subsequent decay into 2 Z bosons 
can produce 4 charged leptons in the final state. This is the dominant background 
process. 

However, using a proper set of cuts one can diminish this background effect, thereby, en- 
hancing the signal. The cuts used in our analysis can be listed as follows: 

• Each of the produced charged lepton (e or // type) must have the pseudorapidity 77 
satisfying — 3 < 77 < 3. 

• Each of the produced charged lepton must have transverse momenta pt greater than 
15 GeV. 

After applying all the above cuts, the cross section for the process pp{gg) ^ ZZ 

(the signal) and pp{gg) ^ h ^ ZZ and pp{gg) — > {EW,non — 

resonant) (the backgrounds) are obtained. In Fig. 3 we have plotted these 

cross sections for the signal and background against ms (GeV). For the signal ms — 
and (0) = 250, 500, 750, 1000 GeV with the SM Higgs mass = 185 GeV. For the 
background ms — nih for the resonating Higgs channel, besides other non-resonating electro 
weak channels. We considered the cases in which Z decay into electrons and muons {I — e, /i) 
with total branching ratios 0.0682. In addition, an average detection efficiency of 80% per 
lepton has been assumed. Below, we make the following observations 

• From Fig. 3, we see that there is a substantial enhancement of the rates (signal) in 
comparison to what is obtained in the Standard Model. The bump at around Ms — 
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350 GeV, resembles the top resonance inside the top loop (the scalar (radion/Higgs) 
production channel via gluon fusion). Such a top resonance, encoded in the function 
\{xt) (see Eq.(16)) and appearing in the loop integral, seems to happen at > 2mt 
which is being manifested in the imaginary part of the X{xt) function (Eq.(17)). 

An estimate of the signal and background event rate (obtained by multiplying the cross 
section with the integrated luminosity) is given in Table 1 and 2. 

10 
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FIG. 3 (color online): Plots showing the total cross section of A l{l = e,/x) in pp 
collisions due to radion resonance as a function of the scalar(S) mass ms, with S is 
and h for the signal and the background, respectively. For the signal, we choose 
the SM Higgs mass ruh = 185 (GeV) and {(f)) = 250 (solid curve), 500 (long-dashed 
curve), 750 (short-dashed curve) and 1000 (dotted curve) (GeV). The lower (dashed- 
dotted) curve corresponds to the SM background. We choose = fi^ = s, although 
the result remains more or less unchanged if one chooses = fi^ = m|. 

• From Fig. 3 and Table 1, it is clear that the total cross section (and hence the event 
rate) of the signal falls with the increase in (0), i.e. near resonance it goes as while 
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away from the resonance it goes as (see also the discussion of Sec. 3). 

• Figure 3 (see also Table 1) shows that, for a given (0), the cross section falls with the 

increase in [see Eqs.(20) and (21)]. 

• The total event (signal) rate is greater than that predicted in ADD and RS scenerio 
involving spin-2 gravitons. As has been mentioned before that this enhancement is 
due to two reasons, namely, (a) possibility of having a radion resonance and (b) the 
enhanced coupling at radion-gluon-gluon vertex due to trace anomaly and top loop. It 
is worthwhile to give some number related to the signal cross section which differs from 
the literatures, say, Chaichian et. al. in the ^ = 0, where ^ is the radion-Higgs mixing 
parameter [5] and as we mentioned earlier that such difference is due to the treament 
of the (f) — t — t coupling required in the traingle top loop evaluation. For example, 
below 2mt threshold, say at around — 200 GeV, we find a — ag — 0.047 pb in 
contrast to ~ 0.01 pb (Chaichian et. al) for the radion vev (0) about 1 TeV. Above 
2mt threshold, with the same (0) value, say at = 370 GeV, we find the signal cross 
section a about 0.06 pb, whereas they found 0.025 pb. 



Table 1: Showing the event (signal) rate as a function of radion mass and vev (0). 



m^{GeV) 


(0) (GeV) 




C ifb-'/yr) 


Events/yr 




250 


408 


100 


40800 


250 


500 


100 


100 


10000 




750 


45 


100 


4500 




1000 


22 


100 


2200 




250 


213 


100 


21300 


500 


500 


53 


100 


5300 




750 


24 


100 


2400 




1000 


14 


100 


1400 
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Table 2: Background (bgd) event rate estimation with the assumption that the dominant 
contribution comes from the SM Higgs (with mh — 185 GeV) resonance due to gluon-gluon 
fusion. 





£ ifb-'/yr) 


Events (bgd) /yr 


4.4 


100 


440 



• A bump at around = 350 GcV corresponds to the top resonance inside the top 
loop. Above this resonance (the 2mf threshold), the cross section is mainly dominated 
by the top loop. Note that, before this resonance, the large signal cross section was 
mainly due to trace anomaly, which remains important even after this 2mt threshold. 

• While maintaining the fact that the signal cross section be greater than that of the SM 
background, one can even allow larger m,^ values corresponding to lighter (0) choices, 
the feature which will be explored in the next section. 

4.2 Significance contours 

With an integrated luminosity of 100 fb~^, the above rate indicate a rather impressive 
prospect of detecting the radion resonance via the pair- produced Z bosons and it's subse- 
quent charged leptonic decays. To gauge the situation, however, one should also remember 
that the backgrounds, even after applying proper set of cuts, are never totally eliminated 
and one can compare the signal with this background. In Fig. 4, we show the contour plots 
in the plane of radion mass and radion vev (0) corresponding to -^i— Rsb) — 5, 10, 
where S and VB are the number of events corresponding to the signal and the background 
with the above luminosity. In calculating the backgrounds, we have taken both statistical 
and systematic effects into account, assuming the systematic uncertainty to be 2% of the 
total background and adding it in quadrature to the computed background itself. For each 
curve in Fig. 4, the region above the curve is allowed. For example, for an intermediate 
(heavy) m^, say about 250 (500) GeV , we find lower bounds on (0) about ~ 2741 (778) GeV 
corresponding to Rsb — 5 and ~ 2076 (~ 556) GeV corresponding to Rsb — 10. These plots 
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corresponds to finding the signal over a large region of the parameter space and suggests the 
possibility of having a heavier radion mass with a moderately lower radion vev (0), although 
the possibility of having an intermediate radion mass is not ruled out. 

10000 p \ \ \ \ 1 




100 I ^ ^ ' ' 1 

100 200 300 400 500 600 

[GeV] 

FIG. 4 (color online). Contour plot in the plane of and {(j)) corresponding to Rsb = 
5 (solid line) and 10 (dashed line). 

5 Summary and Conclusion 

The enhanced coupling of the radion-gluon-gluon due to top loop and trace anomaly offers 
an immense possibility of finding radion at Large Hadron Collider. We have considered the 
correct radion-top-antitop coupling and fix the radion production rate at LHC We have 
considered the channel in which the radion is produced as a resonance in gluon-gluon fusion 
and then decay to a pair of Z bosons, which subsequently decays into 4 charged leptons, 
l{= e,n). The SM background comprising 4 I final sets is found to be smaller that that due 
to the signal (radion resonance) even after applying a proper set of cuts for a wide range of 
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radion mass (m^) values. For m,/, ranging from 180 to 1000 GeV, the signal cross section 
dominates over that due to SM background. Exploiting the signal and background ratio, 
we make the contour plots in m<^ — {(p) plane. The contour plot suggests the possibility of 
having an intermediate to heavier radion of mass about ~ 250 (500) GeV with (0) about 
2076 (556) GeV at the lOcr discovery level. 
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